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1 REASON  FOR  ISSUANCE 

To  document  and  describe  the  use  of  a large  area  spectral 
yield  model  over  the  USSR.  Primary  development  and  testing  done  by 
Pat  Ashburn  of  the  Foreign  Crop  Condition  Assessment  Division, 
Foreign  Agricultural  Service. 

2 COVERAGE 

This  Technical  Memorandum  describes  and  evaluates  the  operational 
application  of  the  AVI  Yield  Model  for  wheat  in  1979  and  1981. 
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PART  1.0  INTRODUCTION 


1 . 1 PURPOSE 

The  purpose  of  this  paper  is  to  report  on  the  development  and  use  of  a 
spectrally  driven  wheat  yield  model  that  is  derived  from  LANDSAT 
spectral  data.  This  paper  is  a summary  of  how  this  model  is  used 
operationally  and  is  not  intended  to  show  detailed  research  on  how  it 
was  developed. 

1.2  SCOPE 

The  scope  of  this  paper  is  limited  to  the  transformation  of  LANDSAT 
multispectral  (MSS)  data  into  a vegetation  index  number  (VIN)  that  is 
used  to  derive  the  yield  of  wheat  for  the  USSR.  Also  in  the  scope  is 
the  aggregation  of  the  vegetation  index  number  to  the  hierarchal 
levels  of  the  USSR  for  the  prime  purpose  of  deriving  a yield  figure 
for  the  country. 

1.3  BACKGROUND 

The  United  States  Department  of  Agriculture  has  a requirement  to 
provide  current  analysis  of  foreign  crop  production.  This  requirement 
led  to  the  establishment  of  a new  Division  within  the  Foreign 
Agricultural  Service  (FAS)  designed  to  take  advantage  of  satellite 
technology.  This  Division,  the  Foreign  Crop  Condition  Assessment 
Division  (FCCAD),  is  currently  located  in  Houston,  Texas.  Its  primary 
function  is  to  monitor,  assess  and  report  on  the  condition  of  crops 
and  agriculturaly  related  events  in  foreign  countries.  To  satisfy 
these  requirements,  the  FCCAD  has  developed  operational  capabilities 
including  a geogridded  reference  data  base,  an  interactive  image 
processing  system,  and  a rapid  image  transformation  system. 

The  development  of  the  vegetation  index  numbers  in  the  USSR  and  how 
they  are  calculated  are  described  in  detail  in  Technical  Memorandum 
Number  11,  "Utilization  of  The  Vegetation  Index  Number  in  the  USSR," 
published  by  the  USDA/FAS  in  Houston,  Texas,  in  1981  . This  paper  is 
intended  to  describe  the  continuation  of  that  study  and  how  the 
vegetation  index  number  has  been  incorporated  into  the  aggregation 
schemes  developed  for  the  geogridded  data  base. 

Results  from  research  on  spectral  reflectance  and  wheat  yield  by 
J.  Kristian  and  F.  J.  Siddoway  at  J^he  Northern  Plains  Soil  and  Water 
Research  Center  in  Sidney,  Montana*",  are  very  encouraging.  Their 
findings  of  very  high  correlation  of  VINs  to  final  yield  of  wheat 
(r  = .98)  help  provide  some  of  the  basic  knowledge  for  the 
development  and  use  of  a large  area  procedure. 
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PART  2.0  THE  VEGETATION  INDEX  MODEL  FOR  WHEAT  YIELD 


2.1  HOW  IT  IS  CALCULATED  AND  USED 

The  Vegetation  Index  Model  for  wheat  yield  in  the  USSR  is  a spectrally 
driven  yield  model  derived  from  LANDSAT  (MSS)  data  (Figure  1).  It  has 
evolved  over  several  years  of  semioperational  use  in  the  USSR  for  both 
spring  and  winter  wheat.  It  has  now  been  transformed  into  a simple 
linear  equation  that  only  requires  two  inputs,  an  average  calculation 
for  the  Ashburn  Vegetation  Index  number  (AVI)  and  the  average  crop 
calendar  stage. 
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2.2  THE  MANUAL  MODEL 
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Solid  lines  represent  curves  for  which  reported  yields  were  available. 
Dashed  lines  represent  curves  for  which  some  yield  information  was 
reported  but  not  enough  to  be  sure  of  its  exact  location.  By  the  fall 
of  1981  enough  data  had  been  collected  to  determine  the  exact  locations 
of  most  of  the  curves  represented  in  this  yield  model.  Also  a certain 
degree  of  confidence  in  the  use  of  the  model  was  obtained.  A decision 
was  made  to  determine  the  linear  algorithm  for  the  green-up  side  of  this 
curve.  The  intent  here  was  to  see  if  a simple  linear  equation  could  be 
developed  and  programmed  into  the  computer  systems  of  the  FCCAD.  If  so, 
this  would  allow  for  a very  rapid  turnaround  of  yield  information  from 
LANDSAT  data.  Limited  results  of  this  effort  are  shown  in  Table  5. 
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FIGURE  1 - YIN  YIELD  MODEL  FOR  WHEAT 
( QUINTALS /HECTARES ) 
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FIGURE  2 - BMTS  and  Phenological  Growth  Stages 
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2.3  THE  MATHEMATICAL  MODEL 


Table  1 shows  the  crop  calendar 
and  the  yields  (YLD)  that  were 
values  were  fed  into  a Minitab 
that  explained  93.7  percent  of 


(CC),  the  vegetation  index  number 
derived  from  the  manual  model.  The 
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the  variation.  This  equation  is: 
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TABLE  1 - THE  STATISTICS  OF  THE  VIN  MODEL  FOR  WHEAT 
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The  algorithm  was  programmed  into  a TI-59  calculator  and  used  oper- 
ationally during  the  calendar  year  1981.  Limited  results  are  shown  in 
Table  5.  All  of  the  USSR  was  not  evaluated  due  to  a lack  of  data  or  a 
lack  of  time  to  develop  the  data. 
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3.1  THE  ORBIT  AND  DATA 

Several  computer  programs  and  procedures  were  developed  to  provide  a 
total  approach  for  calculating  yields  at  the  Oblast  level  for  the 
USSR.  These  start  with  the  VIN  calculation  for  the  individual  picture 
element  of  the  LANDSAT  data,  and  end  with  a weighted  average  for  the 
Oblast.  Several  pitfalls  must  be  recognized  and  corrected  in  the  use 
of  these  data  and  the  approach.  The  following  is  a sequential  step  by 
step  approach  of  the  procedures. 

Data  are  collected  for  every  odd  numbered  orbit  in  the  USSR  during  the 
growing  season.  This  allows  for  the  collection  of  the  same  spot  on 
the  earth  surface  every  18  days,  which  is  a long  time  between  data 
collection,  especially  if  the  area  is  cloud  covered  on  the  overpass 
day  and  no  usuable  data  for  the  crop  are  collected.  Also,  if  the  area 
is  cloud  covered  during  maximum  green,  the  accuracy  of  the  VIN  calcu- 
lation for  yield  is  degraded. 

When  clear  data  are  sampled  during  the  growing  season  for  wheat,  VIN 
masks  are  created  for  the  full  frames.  A VIN  mask  is  a classifi- 
cation map  that  shows  the  location  of  each  wheat  picture  element  in 
the  scene.  They  are  done  by  using  the  IMDACS  VIN  Processor.  This 
processor  does  not  produce  an  area  perfect  class  map,  but  produces  a 
vigor  map.  This  class  map  or  mask  can  often  be  made  area  perfect, 
given  enough  time  and  enough  data. 

3.2  VIGRDM 

When  the  map/masks  for  wheat  are  created,  a program  called  VIGRDM  is 
used  to  divide  the  images  into  their  appropriate  grid  cells  (Table  2). 
The  average  VINs  for  each  of  the  cells  are  calculated  along  with  the 
amount  of  good  data,  the  percent  of  green  pixels,  the  names  of  the 
crops,  and  the  date  of  the  images.  Only  those  grid  cells  that  are 
totally  within  the  data  are  stored.  Consequently,  out  of  a possible 
16  cells  in  a full  frame  of  data,  no  more  than  5 are  commonly  stored. 
Therefore,  the  data  are  only  samples  and  present  all  the  problems  en- 
countered with  sampling.  However,  the  data  calculated  for  each  of  the 
cells  are  good,  and  hopefully  are  representative  of  the  Oblast. 

3.3  VISUMY 

The  program  used  to  produce  an  aggregated  figure  for  the  Oblast  and 
for  higher  geographic  levels  is  called  VISUMY.  A weighted  average 
using  the  number  of  green  pixels  in  each  of  the  grid  cells  and  their 
average  VINs  are  aggregated  to  the  Oblast  level  to  provide  one  VIN  for 
the  Oblast.  A simple  average  of  the  percent  of  green  in  each  grid 
cell  is  taken  to  provide  one  number  for  the  percent  of  green  pixels  in 
the  Oblast.  These  two  numbers  provide  the  VIN  for  use  in  the  yield 
calculation,  and  the  amount  of  wheat  acreage  in  the  Oblast. 
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TABLE  2 - VEGETATIVE  INDEX  NUMBER  GRID  CELL 
FOR  THE  MASKED  IMAGE  (VIGRDM) 


TABLE  2 - VEGETATIVE  INDEX  NUMBER  GRID  CELL  FCR  THE 
HASHn-IMACE-QapKDK). 


•••  PROGRAM  VIGRDM  ***  MASKED  IJ  CELL  V1N  COMPUTATION 


ENTER  3-BYTE  CROP  MASK  FILE  DUAL  I F 1 EK : 

enter  image  file  name  else  blank  tu  uun . 

INPUT  FILES:  HPR=  [70,7]V172021V0.IMH;31S1  I HG=  [ 7 0 , 7 ] V 1 7 2 02 1 V 0 . I MD 7 3 1 5 1 

IMAGE  1D=V172021V0  LATtI  7 24  bl  NUM  CHAN:  4 NUM  P1X:  510  SKIP:  5 ENTS 

NUM  LINES?  510  SKIP?  5 

CROP  MASK  TO  BE  APPLIED:  [ 7 0 , 7 1 V 1 7 2 0 2 1 V 0 . W HS 

FILES  HAVE  BEEN  OPENED. 

1:  29b  301  J:  191  194 

1:  29b  301  J:  191  194 

Mi*  u:  5 Standard  crop  code:  wht 

1 = 299  J=  1 9 2 COHR:  S SUN:  49  \GOUU=100.0  %AHEA=  55.6  A 23.1  L 31.5  K 9.1 

1 = 299  J=  1 9 3 CORK:  S SUN:  49  %gOod=1oo.O  *AREa=  52. b A 25.0  L 32.5  K 10.5 

ENTER  IMAGE  FILE  NAME  else  blank  to  quit. 

I JNPU1  FILES:  HDR=  170,7) V172023V0.lMH;3lbl  1MG=  [7O,7)V172O23V0.lMD;3l51 

ImaCE  ID=V172023V0  DATE:  7 24  HI  NUM  CHAN:  4 NUM  Fix;  510  SKIP:  5 ERTS 

KUR'LII.ES:  510  SKIP:  5 

CROP  MASK  TO  BE  APPLIED: [ 7 0 , 7 J V 1 7 2 0 2 3 V 0 , WHS 

t 1 LES  HAVE  BEEN  OPENED. 

1:  303  306  J:  1 B 7 190 

1:  303  306  J:  FT?  T90  ' — ' 

NUM  1J:  3 STANDAHD  CROP  CODE:  WHT 

1=304  jiiee  cohr:  s sun:  so  %good=  97.5  \area=  7a.e  a 29.0  l 3o.3  k lb.b 

1 =305  J=  1 8 8 CORR:  S SUN;  50  AG0Ut>=  99.3  *AKEA=  45.2  A 18.5  L 21.7  K 9.b 

1=305  J=189  cork:  S sun:  50  *GOOD=  99.5  AAREA=  27.6  A 15.5  L 18.9  K 7.3 
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One  product  from  this  program  is  shown  in  Table  3.  In  this  table, 
special  attention  should  be  paid  to  the  number  of  grid  cells  used  to 
calculate  the  average  VIN  for  each  level  of  the  geographic  hierarchy 
and  the  date  for  each  calculation.  Also  note  that  the  winter  grains 
masks  (MGW)  are  separate  calculations  of  the  data.  These  masks, 
whether  for  winter  grain,  rye,  or  wheat,  carry  the  numbers  used  in  the 
VIN  wheat  yield  model.  Also  note  that  only  a few  numbers  are 
calculated  for  specific  crops.  Few  masks  were  produced  because  of 
time,  poor  data,  or  the  date  of  the  data.  The  data  elements  in  Table 
3 are  : 


1 . Calendar  date 

2.  18  day  period 

3.  Number  of  grid  cells  used  in  the  calculation 

4.  Satellite  type 

5.  Type  of  crop  mask 
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not  difficult  when  all  the  wheat  in  a given  Oblast  is  in  a similar 
condition,  but  in  large  Oblasts  with  varying  crop  conditions,  this 
becomes  more  difficult. 
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yield  of  wheat  for  e 


to  determine  the  accuracy  of  aggregating  the 
Is  to  higher  levels.  This  could  be  done  but  is 
this  paper.  The  average  VINs  for  the  aggregated 
the  calculated  crop  calendar  to  determine  the 
ach  area.  These  are  summed  to  the  country  level. 
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TABLE  3 - VEGETATIVE  INDEX  NUMBER 
SUMMARY  (VISUMY) 
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PART  4.0  RESULTS 

4.1  PAST  RESULTS 

Past  use  of  the  manual  model  has  provided  accurate 
These  results  were  derived  from  data  taken  over  the 
of  the  New  Lands  in  1979.  The  information  derived 
was  used  to  provide  some  additions  and  revisions  to 
model.  This  allowed  for  a better  understanding  of 
between  green  biomass  and  yield.  Also,  relationshi 
and  conditions  existing  after  the  index  numbers  wer 
better  understood. 

It  became  very  evident  that  additional  systems  were 
mine  yield  detracting  characteristics  after  the  veg 
numbers  were  last  calculated.  For  instance,  in  197 
was  directly  attributed  to  a late  harvest.  This  co 
accuracy  determination  of  the  model,  since  the  mode 
late  harvest  had  not  been  independently  tested. 

4.2  CURRENT  RESULTS 

In  Table  5,  the  Ukraine  SSR,  the  North  Caucasus  Region,  and  Kazakhstan 
SSR  are  shown  with  their  average  VINs.  The  estimated  yield  derived 
from  the  VIN,  the  average  VIN,  and  the  reported  yield  from  the  USSR 
are  shown.  It  is  evident,  from  looking  at  this  table,  that  little 
yield  information  is  available  for  this  past  year.  This  hampers  an 
accuracy  determination  for  the  model.  Consequently,  more  time  is 
required  before  enough  yield  data  for  a good  evaluation  of  the  model 
and  the  aggregation  procedures  are  completed. 
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TABLE  4 - ESTIMATED  SPRING  WHEAT  YIELDS 
DERIVED  FROM  VINS  FOR  1979 


OBLASTS 

1970-1974 

AVG.  YIELD 

ESTIMATED 
YIELD  SPRING 
WHEAT  1979 

REDUCTION 

DUE  TO  LATE 

HARVEST 

*REPORTED  YIELD 

IN 

TOTAL  GRAIN 

Altay 

13.3 

13.1 

13.4 

Bashkir 

14.1 

19.3 

17% 

16.0 

15.9 

Kemerovo 

13.8 

18.0 

5% 

17.6 

Kokechetav 

12.1 

(11. 4)1 

18.7 

6% 

17.6 

14.9 

Kurgan 

16.0 

( 6 . 6 ) 1 

19.6 

11% 

17.6 

15.0 

Kustanay 

9.7 

(8. 5)1 

17.3 

16% 

14.5 

14.0 

Novosibirsk 

13.0 

19.3 

6% 

18.3 

Omsk 

14.5 

19.3 

11% 

17.3 

17.9 

Orenburg 

11.5 

11.0 

13.5 

Pavlodar 

7.7 

( 6 . 6 ) 1 

11.0 

12.2 

Severo-Kazakh 

14.6 

( 1 2 . 3 ) 1 

20.2 

10% 

18.0 

18.3 

Tselinograd 

10.0 

(9 . 9 ) 1 

15.5 

17.0 

*MOST  OF  THE  CROPLAND  IS  SOWN  TO  SPRING  WHEAT 
1 TOTAL  GRAIN  1971-1975 
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TABLE  5 - VIN  NUMBER  YIELD  AND  RESULTS 


WINTER 

WHEAT 

GEOGRAPHIC 

EST.  YIELD 

REPORTED 

AREA 

AVERAGE 

VIN 

FROM  VIN 

YIELD 

Western  Ukraine 

37.5 

27.1 

Donets-Dnepr 

40.9 

28.5 

Southern  Ukraine 

29.0 

21.8 

Zaporozhye  Oblast 

47.5 

32.3 

32.7 

Moldavia  SSR 

29.7 

21.9 

Kransnodar  Kray 

48.6 

34.2 

Stavropol  Kray 

46.8 

31.8 

32.8 

Rostov  Oblast 

35.3 

25.2 

23.0 

N.  Caucasus  Region 

42.1 

28.9 

SPRING 

WHEAT 

Uralsk  Oblast 

13.6 

9.3 

Turgay  Oblast 

20.4 

15.0 

Tselinograd  Oblast 

21.4 

15.6 

Kokchetav  Oblast 

20.3 

16.0 

Severokazakh  Oblast 

18.8 

16.2 

Pavlodar  Oblast 

7.5 

7.0 

Karaganda  Oblast 

16.6 

12.0 

Semipalat insk  Oblast 

8.1 

9.3 

Kazakh  SSR 

18.9 

13.2 

mi: 

O O L. 
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CONCLUSION 


1 CONCLUSIONS 

This  is  the  first  time  that  the  VIN  wheat  yield  model  has  been  tested 
over  such  a large  area  of  the  USSR.  And  as  of  this  writing,  the 
results  are  inconclusive,  mostly  due  to  a lack  of  reported  yield  data. 
This  does  not,  however,  alter  the  past  successful  results  of  the  use 
of  the  model  over  smaller  areas  of  the  USSR,  as  shown  in  the  data  from 
1979.  Sufficient  results  were  obtained  during  1979  and  the  past  two 
years  to  believe  that  the  model  could  be  put  into  use  on  a wide  scale. 
This  of  course  requires  that  data  be  available,  and  that  a sufficient 
number  of  resources  be  put  into  the  classification  and  calculation  of 
the  appropriate  VIN  data  for  countries  of  interest.  It  seems  appro- 
priate for  AgRISTARS  to  include  the  United  States  in  a test  project, 
so  that  more  data  for  county  or  state  size  areas  would  be  available. 
This  activity  would  also  help  determine  the  exact  location  of 
additional  curves  within  the  model. 

The  VIN  model  can  provide  a one  time  capability  for  deriving  the  yield 
over  a very  large  area,  the  best  time  being  at  or  just  before  maximum 
greenness  of  the  plant.  This  accounts  for  the  accumulation  of  yield 
up  to  that  image  date  but  cannot  account  for  any  yield  affecting 
events  after  the  date.  So  additional  data  (weather)  should  be  used  to 
monitor  any  adverse  effects  after  the  last  date  of  imagery.  It  is 
also  apparent  that  the  second  component  of  this  vegetation  index 
number  yield  model,  the  crop  calendar,  be  as  accurate  as  possible. 


Only  time  will  tell  how  accurate  the  procedure  and  the  model  are  for 
estimating  wheat  over  large  areas  of  the  USSR.  This  paper  has  not 
addressed  the  accuracy  problem,  but  has  provided  a model  and  an 
aggregation  procedure  with  existing  software  that  could  be  used  as  an 
analytical  tool  in  other  countries  as  well  as  the  USSR. 
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GLOSSARY  AND  REFERENCES 


AVI  ~ Ashburn  Vegetation  Index  Number 

IMDACS  - Integrated  Multivariate  Data  Analysis  and  Classification  System 

KVI  - Kauth  Vegetation  Index  Number 

LAI  - Leaf  Area  Index  Vegetation  Index  Number 

MASK  - A Classification  Map  of  Nonthresholded  Picture  Elements 
MSS  ~ Multispectral  Scanner  from  the  LANDSAT  Satellite 
VIGRDM  - Vegetation  Index  Grid  Cell  Map  Program 

VIN  - Vegetation  Index  Number  - A numerical  number  that  describes  the 
greenness  of  the  picture  element 

VISUMY  - Vegetation  Index  Number  Summary  Program 


REFERENCES 


1.  Ashburn,  Pat,  "Utilization  of  the  Vegetation  Index  Number  in  the  USSR," 

TM-11  (May  18,  1981),  The  Foreign  Crop  Condition  Assessment 
Division,  FAS/USDA. 

2.  Aase , J.  K. , and  Siddoway,  F.  H. , "Spring  Wheat  Yield  Estimates  from 

Spectral  Reflectance  Measurements,"  Northern  Plains  Soil  and  Water 
Research  Center,  Sidney,  Montana,  January  1,  1981. 


El-1 


9 0 


1982 


TM-14 


EXHIBIT  2 
CPARA,  2.3) 


TABLE  1 - THE  STATISTICS  OF  THE  VIN 
MODEL  FOR  WHEAT  -- 


— tLi — 

THE  REGRESSION  EQUATION  IS 

Y = 

20.0  - 5 

.01  XI  ♦ 0.545  X 2 

COLUMN 

COEKFICIL'NT 

ST.  DEV. 

OF  COEF.  ’ 

r^rrsnu-^ 

COEF/S.D. 

XI 

CC 

20.0175 

-5.0095 

0.4510 

0.1932 

44.38 

-25.92 

n 

~TVJ 

"8T5  4 55 

0.D095 

57T27 

I HE  SIT 

s = 

~de v . ur~ 

1.430 

V ABOUT  regression 

LINE  iS 

STTfTT- 

Ti6-  ' 3 T~= 

233 LEG REE 5 OF  FREEDOM 

H-SUUAHED  = 93.7 

R-SQUARED  = 93.6 

PERCENT 

PERCENT,  ADJUSTED 

FOR  D.F. 

ANALYSIS  OF  VARIANCE 


DUE  TO 

DF 

SS 

MS 

=SS/DF 

RLCReSSIOnI  7 

RESIDUAL  233 

7053.319 

476.506 

3526.090 

2.045 

TOTAL 

7T5 

'7529.964 

FURTHER 

ANALYSIS 

OF  variance 

SS  EXPLAINED  BY 

EACH  VARIABLE  WHEN 

ENTERED  IN 

the  order  given 

DUE  TO  DF 

REGRESSION  2 

SS 

7053.379 

CC 

AVI 

1 

1 

345.149 

6708.230 

XI 

T 

PRED.  r 

ST. DEV. 

ROW 

97 

CC 

2.80 

YLO 

8.0000 

VALUE 

10.9002 

PRED.  Y 
0.1661 

RESIDUAL 

-2.900? 

ST. RES. 

-2 . 04W 

146 

147 

3.00 

3,40 

21.0000 

8.0000 

24.0805 
8 .4400 

0.1292 

0.2372 

-3.0805 

-0.4400 

-2.  1 6 R 
-0.31  X 

140 

109 

3.40 

3.40 

8.5000 

33.0000 

8.9854 
J1 .3497 

0.2292 

0.2359 

-0.4854 
t .6503 

-0.34  X 
1.17  X 

190 

.191 

3.40 

3.40 

33.5000 

35.0000 

31.8951 

32.440b 

0.2440 

0.2522 

1.6049 

2.5594 

1.14  X 

1.82  X 

192 

.193 

3.40 

3 .40 

35.8000 

36.8000 

32.9861 

33.5316 

0.2604 

0.2688 

2.8139 

3.2604 

2. OORX 

2 . 33RX 

194 

195 

3.40 

3 .60 

37.5000 

1 1 .0000 

34.0770 

9.6199 

0.2773 

0.2387 

3.4230 

1.3801 

2 . 4 48  X 

0.98  X 
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233 

3.60 

3.60 

1 1.7000 
32.0000 

10.1654 

30.3478 

0.2313 

0.2293 

1.5346 

1.6522 

1.09  X 

1.17  X 

234 

235 

3 .60 
3.60 

33.0000 

34.0000 

30.8932 

31.4387 

0.23fo6 

0.2441 

2.1068 

2.5613 

1.49  X 

1.82  X 

236 

3.60 

35.0000 

31.9842 

0.2517 

3.0158 

2. 1 4 R X 

R DENUTES  AN  OBS'.”wiTH  A LARGE  ST,  RF.S. 

X DENOTES  AN  UBS.  WHOSE  X VALUE  GIVES  IT  LARGE  INFLUENCE. 
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--  PkllJT  Cl  C2  C3  C 5 

COLUMN 

CC 
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YLD 

ESTYLD 

COUM 

216 
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2T6 

2TS 

HOw 

1 
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3. 

1 0 .0000 
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4. 
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4 
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1 

1 .600013 
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1 .80000 
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1 7 . 0000 

15. 3b42 

1 
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~ 57 
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e 

1 .80000 
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1 9 . 0000 

16.4551 

y 

1 . E C 0315 

u. 

1 9.  50TJ73 
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JO 

1.80000 

12. 
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17.5460 

ii 

1.80000 

13. 

20.0000 

18.0915 

12 

1 .80000 

14. 

20.0000 

18.6370 

13 

1 . 8 0 0 0^ 

15. 
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1 9 . 1775 

14 

1 .80000 

16. 
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13 

1 . 80000 
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20.8OOO 
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16 
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18. 
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n 
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n 
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21 
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2.00000 
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19.9000 

18.4477 

9? 
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22. 
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93 
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23. 
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94 
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24, 
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20.0000 

21.1750 

97 

2.80000 

9. 

8 . OOOO 

10.9002 

98 

2 .80000 

10. 

9 . 0000 

11.4456 

99 

2.80000 

1 1 . 

10.0000 

11.9911 

100 

2 .80000 

12. 

11. 0000 

12.536b 

101 

2.80000 

13. 

12.0000 

13.0820 

1 02 

2.80000 

14. 

14.0000 

13.6275 

103 

2 .80000 

15. 

15. 0000 

14.1730 

1 04 

2.80000 

16. 

16.0000 

14.7185 

1 05 

2 .80000 

17. 

lb. 3000 

1 5 . 2b  39 

106 

2 .80000 

18. 

1 6 „ bOOO 

15.8094 

107 

2.80000 

19. 

16.8000 

16.3549 

1 OR 

2.80000 

2_0j 

1 7 . 0000 

16.9003 

1 09 

2 .80000 

21. 

1 7 . 3000 

17.4458 

no 

2 .80000 

22. 

1 7 . bOOO 

17.9913 

in 

2 .80000 

23. 

17 .8000 

18.5367 

1 12 

2.80000 

24. 

1 8 . 0000 

19.0622 

1 1 3 

2.80000 

25. 

1 9 . 3000 

19.6277 

1 1 4 

2.80000 

26, 

19.6000 

20.1731 

115 

2 .80000 

27. 

19.8000 

20. 7 1 b6 

1 1 6 

2 .80000 

2b. 

20 . 0000 

21.264] 

1 1 7 

2.80000 

29. 

20.2000 

21 .8096 

1 1 P 

2.80000 

_3Q_. 

20, 5000 

22.3550 

1 1 9 

2.80000 

31 . 

2 1 .0000 

22.9005 

L?o 

3.00000 

9 . 

7.5000 

9.8983 

121 

3 .00000 

10. 

8.0000 

10.4437 

1 22 

3 .00000 

1 1 . 

8.5000 

10.9892 
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EXHIBIT  2 
(PARA,  2.3) 


--  PkIIJT^CTc2' 

COLUMN  CC 

AVI 

YLD 

E5T  Y LD 

C u l)  M 

HUM 

236 

236 

23"5 

236 

123 

3.00000 

12. 

9.0000 

11.5347 

i 124 

3 . 00000 

13. 

1 0 . 0000 

12.0601 

125 

3T7H5000 

14. 

12.0055 

12.625b 

1 26 

3.00000 

15. 

1470000 

13.1711 

127 

3.00000 

16. 

T5. 000(5 

T3  .1 1 66 

126 

3.00000 

17. 

15.5000 

14.2620 

~TT9~ 

3.00000 

16. 

1 6 . &OO0 

147  8 07  5 

130 

3.00000 

19. 

16.3000 

15.3530 

13l 

3.00000 

20. 

1 6 . 6000 

1 5 . 8 9~54 

132 

3.00000 

21. 

1 6 .9000 

16.4439 

133 

3.00000 

22. 

1 7 . 0000 

16.9894 

134 

3.00000 

23. 

17.2000 

17.5348 

135 

3.00000 

“24T“ 

17.5000 

16.0603 

136 

3.00000 

25. 

17.7000 

18.6258 

137 

3.O0OOO 

26. 

17.9000 

19.1713 

136 

3.00000 

27. 

18.0000 

19.7167 

139 

3.00000 

28. 

18.5000 

20.2622 

140 

3.00000 

29. 

19.0000 

20.8077 

. m 

3 .00006 

30. 

l y. booo 

— 21.3531 

142 

3.00000 

31. 

20.0000 

21.8986 

143 

3.00000 

32. 

20.3000 

22.4441 

144 

3.00000 

33. 

20.6000 

22.9895 

145 

3 .ooooo 

34. 

20.8000 

23.5350 

146 

3 . 00000 

35. 

21 .0000 

24 . 0805 

147 

3.40000 

757 

6.0000 

b . 4450 

146 

3 .40000 

1 1 . 

8 . 5000 

8.9654 

149 

3.4  00IT5 

12. 

977)0(15 

9.5309 

1 50 

3.40000 

13. 

10.0000 

10.0764 

151 

J7T00CT0 

141 

— I r.  D 0 0 0 

I 0762  I S 

152 

3 .40000 

IS. 

12.0000 

1 1 . Ib73 

T51 — 

3 . 4 0 0 0 0 

16. 

13.0000 

11.7126 

1 64 

3.40000 

17. 

1 4 . 0000 

12.2582 

155 

J74  DOTTC 

18. 

15.0000 

12. 8 0 3 7 

156 

3 .40000 

19. 

15.5000 

13.3492 

1 57 

3T4  000TS 

20. 

1 5 .6000 

13.694b 

1 66 

3 .40000 

21  . 

1 6 . 0000 

14.4401 

1 5y 

3.4G0&0 

'"'227~' 

16.3000 

14.965b 

160 

3.40000 

23. 

16.6000 

15.531 1 

161 

3.40000 

“ 24  T 

' J 77T0  0 <35 

16.0765  ” 

162 

3.40000 

25. 

17 . 3000 

16.6220 

163 

37  40060  ' 

7b~. 

17.5000 

17.1675 

164 

3.40000 

27. 

17.7000 

17.7129 

1 65 

“ 3 .40000 

2b. 

17. 9000 

“1672564 

166 

3.40000 

29. 

1 8 .0000 

18.8039 

167 

3.4  0TTtrO 

"30. 

“I  6775055 

19.3493 

168 

3.40000 

31. 

19.0000 

19.8948 

16  9 

3740000" 

32. 

975000” 

2 074  403 

170 

3.40000 

33. 

20 . 0000 

2 0 . 9 b 5 8 

171 

3 .40000 

347 

20.5000 

21.5312 

172 

3.40000 

35. 

2 1 .0000 

22.0767 

173 

3.40000 

567 

21.5000 

2 27  62 22  ' 

174 

3.40000 

37. 

22 . 0000 

23.167b 

175 

3.40000 

38. 

2279000 

23.7137“ 

176 

3.40000 

39. 

23.5000 

24.2566 

177 

3.40000 

40. 

24.1005 

2 4.H04  0 

178 

3.40000 

41  . 

24 . 8000 

25.3495 

' 179 

3.40OOO 

477 

— 25 . 300o 

"25.6950 

180 

3 .40000 

43. 

25.9000 

2b. 4404 

1 b 1 

3 .40000 

44  . 

26.5000 

26.9659 

162 

3.40000 

45. 

27 . 0000 

27.5314 

163 

3 . 4 0 0 075 

4b. 

“26 .0000 

“2  67  0769 

1 64 

3 . 4 0000 

47. 

28.6000 

26.6223 

1 b5 

3 . 40000 

4b. 

29.2000 

29. Ib76 

l b6 

3. 40000 

49. 

30.0000 

29.7133 

167 

3 . 40000 

50. 

30.5500 

30.2587 

1 86 

3.40000 

51  . 

32 . 5000 

30.6042 

E2  4 
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EXHIBIT  l 
(PARA.  2.3) 


ACftICULTURAL  library 


1022378739 


— =-  pk  niTTTTT- c 

COLUMN  CC  AVI 

Y LD 

ESTYLD 

CORK! 

KOM 

236 

236 

2 31 

236 

IBS)  . 
1 90 

3 .40000 

3 .40000 

52. 

5 3^. 

33.0000 

33.5000 

31.3497 

31.8951 

191 

192 

3.40000 

3.40000 

54. 

_ _ 55. 

56. 

57. 

35.0000 

3S.H000 

32.440b 

32.9861 

193 

194 

3.40000 

3 .40000 

36.0000 

37.5000 

33.531b 

34.0770 

195 

3.60000 

14. 

11. 0000 

9.6199 

1 9fc 

3.60000 

15. 

1 1 .7000 

10. 1654 

197 

3.60000 

16. 

12.0000 

10.7109 

196 

3 .60000 

17. 

1 4 , 0000 

11.2563 

1 99 

3.60000 

IB. 

1 4 . 5000 

1 1 .0016 

200 

3.60000 

19. 

1 5.0000 

12.3473 

201 

3.60000 

20. 

15.3000 

12.0928 

207 

3.6QPQO 

21. 

15.6000 

13.4382 

203 

3.60000 

22. 

16.0000 

13.9037 

204 

3 .60000 

23. 

16.3000 

14.5292 

205 

3.60000 

24. 

16.6000 

15.0746 

206 

3 .60000 

25. 

17.0000 

15.6201 

207 

3.60000 

26. 

17.2000 

16. lbbfa 

706 

3.60000 

27. 

17.5000 

16.7110 

209 

3.60000 

2 B . 

17.7000 

17.2565 

210 

3 .60000 

29. 

1 8 .0000 

17,8020 

211 

3.60000 

30. 

1 fa . 3000 

18.3475 

212 

3 .60000 

31  . 

1 b . 6000 

18.6929 

213 

3.60000 

32. 

19.0000 

19.4364 

214 

3.60000 

33. 

19.5000 

19.9639 

215 

3.60000 

34. 

20.0000 

20.5293 

216 

3.60000 

35. 

20.3000 

21.0748 

717 

3.60000 

3b. 

20.6000 

21  . b 2 0 3 

2 1 6 

3.60000 

37. 

20.9000 

22. 1657 

219 

3.60000 

38. 

2 1 . 0000 

22.7112 

220 

3.60000 

39. 

21 .5000 

23.2567 

221 

3.60000 

40. 

22.0000 

23.8021 

222 

3.60000 

41. 

22.5000 

24.3476 

223 

3.60000 

42. 

23.0000 

24.8931 

• 224 

3 .60000 

*3. 

23.5000 

25.4386 

225 

3.60000 

44. 

24 . 0000 

25.9o40 

226 

3 .60000 

45. 

25.0000 

26.5295 

227 

3.60000 

46. 

2b. 0000 

27.0750 

226 

3 . 60000 

47. 

27 . 0000 

27.6204 

229 

3.60000 

4 0. 

2«  . 0000 

28.1659 

230 

3 .60000 

4 9 , 

29 . 0000 

28.7114 

231 

3.60000 

50. 

30.0000 

29.2568 

237 

3 .60000 

51. 

3 1 . 0000 

29.8023 

233 

3.60000 

52. 

32.0000 

30.3478 

234 

3 .60000 

53. 

33.0000 

30.8932 

235 

3.60000 

54  . 

34 . 0000 

31.4387 

236 

3.60000 

55. 

35.0000 

31  .9642 

nasa-jsc 
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